Failure of the local density approximation in time-dependent 
spin density functional theory 
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It has been known for some time that the exchange-correlation potential in time-dependent density 
functional theory is an intrinsically nonlocal functional of the density as soon as one goes beyond the 
adiabatic approximation. In this paper we show that a much more severe nonlocality problem, with 
a completely different physical origin, plagues the exchange-correlation potentials in time-dependent 
spin-density functional theory. We show how the use of the spin current density as a basic variable 
solves this problem, and we provide an explicit local expression for the exchange-correlation fields 
as functionals of the spin currents. 
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For many years the local density approximation (LDA) 
has provided the much needed handle on the difficult 
problem of approximating the density dependence of the 
exchange-correlation (xc) potential - the single particle 
potential that incorporates the many-body effects in the 
Kohn-Sham equation for the ground state density Q. In 
LDA, the xc potential V xc (f) is simply a function of the 
local density n(r). This approximation is not unreason- 
able as long as the functional derivative of V xc (f) with 
respect to n(f) - the so called exchange- correlation ker- 
nel / IC (f, f) = S Jn{^ - is a sufficiently short-ranged 
function of the distance \ f — f > \ 0. 

However, much recent work |, ^, ^, |7| has demon- 
strated that the requirement of short-rangedness is not 
always fulfilled in physical systems, and when this hap- 
pens the local density approximation is flawed. This does 
not mean that a local description of exchange and correla- 
tion is absolutely impossible, only that such a description 
cannot be achieved in terms of the particle density. 

For example, in the density-functional theory of crys- 
talline insulators, it has been found ^, [| || that the xc 
potential has an "ultranonlocal" dependence on the den- 
sity, due to the fact that the Fourier transform of the xc 
kernel f xc (k,k) diverges as 1/fc 2 for k — > in these sys- 
tems. But, the ultranonlocality disappears if one refor- 
mulates the theory in terms of the electric polarization 
P(r) and the exchange-correlation electric field E xc {r) 
associated with it. 

A similar phenomenon was discovered in the time- 
dependent density functional theory (TDDFT) (|] follow- 
ing the realization that the frequency-dependent LDA^] 
fails to satisfy Kohn's theoremjn], f[j]| . The pathology 
was traced to a singularity of the form in the xc 

kernel f xc (k,k',u) for k —* at finite k' and u. The 
ensuing nonlocality problem was solved by upgrading to 
time-dependent current-density functional theory (TD- 
CDFT), where the basic variable is the current density, 
and its conjugate field is a vector potential@. TDCDFT 
has since been applied to the calculation of the optical 
spectra of solids |TJ] and the polarizability of long poly- 
mer chains [Ol with considerable success. 



In this Letter we show that the nonlocality problem 
occurs in an aggravated form in the time-dependent spin 
density functional theory or, more generally, in the time- 
dependent DFT of any multi-component system. The 
novel features of the spin-dependent problem stem from 
the fact that the xc kernel presents a divergence even 
in the homogeneous electron liquid. More precisely, it 
can be shown that the Fourier transform of the spin- 
dependent exchange-correlation kernel f xc ,aa , {i' — r' ,t — 
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long- wavelength expansion 



in a homogeneous electron liquid has the 
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+ B aa ,(uj) + 0(k 2 ) , (1) 



where A(cu) and B aa >(uj) are complex functions of fre- 
quency, n a is the density of c-spin electrons [a = +1 for 
|-spin and a = — 1 for J,-spin), and n — n-\ + is the 
total density. Since the xc potential created by a small 
density variation 8n a {k,uj) is given by the formula 



v xc ^{k,uj) = y^j xc 



(2) 



we see that Eq. (|l|) rules out the possibility of a local 
connection between V xc . a {f,t) and 5n rT '(r' ,t'). 

The existence of the long-wavelength singularity in 
fxc.a-cr'(k,uj) has been known for some time. It was first 
pointed out by Goodman and Sjolander Q| that the 
third moment sum rule for the spin-density response 
function implies such a singularity. Approximate formu- 
lae for f xc _(k,uj) = f X c,rt(k,w) - f xc ji(k,uA exhibiting 



5[ and |]16[. More re- 
lave shown that, at 



the singularity were proposed in 
cently, D'Amico and Vignale (|l 
low frequency and finite temperature, the singularity is 
related to the friction that arises between up- and down- 
spin currents when they have different average velocities 
(the so-called spin-drag effect). 

By contrast, the implications of Eq. ([!]) for spin den- 
sity functional theory have not been explored so far. This 
is understandable. The singularity (p arises only at fi- 
nite frequency (A(Q) = 0) and therefore does not af- 
fect the static spin DFT. Furthermore, the singularity 
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does not show up as long as one is interested only in 
the density response of spin-compensated systems, since, 
in that case, the relevant combination of xc kernels is 
Zw n <r n a' fxc.aa', which is non-singular. It is only in 
the time-dependent spin DFT JL8J that the issue of the 
singularity becomes really critical not only to the calcu- 
lation of the spin response, but even to the calculation of 
just the density response Q. 

In this Letter we propose a resolution of the nonlo- 
cality problem based on the use of the spin components 
of the current density j j (r , uf) and j j (r , uj) as basic vari- 
ables. We provide an explicit expression for the spin- 
dependent exchange-correlation field E xc , a (r,uj) as a lo- 
cal linear functional of the currents j a . 

The general method for upgrading from the density 
to the current-density formulation is described in de- 
tail in Ref. ]2Q| , so we mention only the essential steps 
here. We introduce a spin-dependent xc vector poten- 
tial A XCtCr (k,oj) (whose time-derivative, iwA xc ^ a {k, uj) = 
E xc (k,uj), is the xc electric field), and notice that this is 
linearly related to the currents in the following manner 

AZ^fcu) = (3) 
uj — 

u' 

where the superscript a denotes the longitudinal (a = L) 
or transverse (a = T) component of a vector relative to 
the direction of k. It is not difficult to see that the lon- 
gitudinal xc kernel defined in this manner coincides with 
the xc kernel introduced in Eq. (|IJ) . The extra factor ^ 



in Eq. (||) exactly cancels the small-A; singularity of f xc , 
and leads to a theory that admits a local approximation. 
The imaginary part of the current xc kernel f xc _ aa i(k, of) 
is expressed in terms of a causal response function as 
follows: 

SmfZ^Abu) = v 1 , 2 am((F a "(fc);F», (-£)))„ , 
vn a n a ik' L 

(4) 

where ({A;B)) U = -j? '{[A(t), B])e lult dt is the linear 
response function associated with the operators A and 
B; = -S%-[H,j$(%)] is the time derivative of the 

Fourier transform of the current-density operator j a (k), 
H is the Hamiltonian, and V is the volume. 

Once the imaginary part of f xc aa .,(k,oj) is known, its 
real part is determined by the Kramers-Kronig dispersion 
relation 

Kc/^Cfe.w) = f^A^oo) 

2^ r, yzmf^Ak,"') 



where V denotes the principal part integral, and the infi- 
nite frequency limit of f xc ,oo' is determined by the third 
moment sum rule. In a three-dimensional electron liquid, 
this sum rule gives 



/;. ^ ' * • ^) ^° [.9Tl(0) - 1] oa> + ±- [a a t ca 8 ary/ + b a e pot , aa ,\ (6) 



where a L = 2, a T = 2/3, b L = 4/15, and b T = -2/15. 
Here <?jj,(0) is the pair correlation function for antipar- 
allel spin electrons at zero separation, t ca is the average 
correlation kinetic energy of the cr-spin component, and 
t V ot,oo' = J dr^[g aa >(r) - 1] is the potential energy 
associated with the interaction between a- and cr'-spin 
electrons. Note that the result for the longitudinal case 
was first obtained in Rcf . [[14j . 



It is evident from the above equations that both the 
longitudinal and the transverse kernels exhibit A- singu- 
larities, which are "cured" by the ^ factor of Eq. (||) . In 
particular, substituting the small-fc expansion F"(k) — 
F?(0) + 0(k) in Eq. (@), where F?(0) is the operator of 
the total force acting on tr-spin electrons, and noting that 
terms of first order in k vanish by inversion symmetry, 



we see that the xc-kernels have the small-fc expansion 
fc^o A(oj) aa'n 2 



fxc,aa'( k ^) 

where 



k 2 An a n 



+ B?A">) + 0(k*), (7) 



and 



$smA(u)) 



^ieA(uj) 



Vn 2 



$Sm((F?;F?)) u 



(8) 



Wire 2 



[<m(o)-i] 



ly r^ mA ^. (9 ) 



uj 2 - uj' 2 



The factor aa' in Eq. (Q) arises from the fact that the 
total force F^ + F^ vanishes by translational invariance, 
so that ((F CT ;F CT /)) W = -aa'{{F r ;F L )) u . Notice also that 
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A{ui) is independent of the direction a - longitudinal or 
transverse. The microscopic expression for B" , is more 
complicated: a simple approximation for this quantity 
will be presented below. 

Substituting the expansion (@) in Eq. (j^) , calculations 
similar to those described in [g0| lead us to the following 
local approximation for the xc field in terms of the spin 
currents 



eE xc , a (uj) = -VV x L c D a A + — V- SwM 

n n 



in 2 A{u>) 



4u> 



(10) 



Here the r dependence has been left implicit, and the xc 
stress tensor cr xc (u)), as well as A(w), is a function of the 
local spin densities, as discussed below. 

Eq. ( p0| ) is the central result of this paper. The first 
two terms on the right are well known: they are, re- 
spectively, the adiabatic LDA contribution and the visco- 
elastic force term, where the stress tensor a xc . a (uj) is re- 
lated to B xc arT i by obvious extensions of the formulae 
reported in [pof . The expression for the xc stress tensor 
is 



( du a 



V drj 



= £ 

a' 

+ Cxc,<tct'V ■ U a i5; 



dn 



where u a = j a /n a i and 

Vxc,aa' = ; B aa ,{uj) , 



(11) 
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FIG. 1: Imaginary part of A(w) evaluated from Eq. (^) with 
the correction factor given in Eq. ([Hi]). The values of a(r s , 0) 
are 1.92, 3.36, and 7.49 at r s = 1,2, and 4 respectively. 
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FIG. 2: Real part of A(ui) calculated from Eq. 



Cxc,™< = -^fcH - jb^M - 4' c ,^] , (is) 

where e" „„, = -J^r^- The last term in Eq. (pi) 

is new, and comes directly from the p- singularity of 
Eq. (Q). The essential feature of the new term is that 
it produces damping of the spin current proportional to 
the relative velocity between up- and down-spin elec- 
trons. This makes it readily distinguishable from the 
usual viscous friction contained in the second term, which 
is proportional to the derivatives of the velocity field. 
The physical reason for the difference is that, whenever 
up and down spin currents travel with different average 
velocities, they exert friction on each other: the "spin 
drag coefficient" is 7(cj) = , m A ^ . Of course, like 

all the quantities considered here, j(u>) is complex and 
frequency-dependent, and, in the limit of zero frequency, 
its real part can be shown to be related to the spin diffu- 
sion constant D s by the Einstein relation D s — ^^ttvt, 

where \s is the static, macroscopic spin susceptibility. 

Unfortunately, an exact calculation of A(u>) from 
the microscopic expressions (^) and (^) is beyond the 



reach of present-day many-body techniques. However, 
we can obtain a rather good approximation with the 
help of the following exact results: (i) For co — ► 0, 
^smA(uj) oc u! 3 and 5ReA(w) cx ui 2 ; (ii) For large u>, 
ZmA(u;) -> and ^ eA{u)) 

_l$ f ?mni [gn{0) _ lL Here Q = where Ep] 

is the Fermi energy for majority spin electrons and 
£ = " 1 n " 1 measures the degree of spin polarization, 
and a = (4/97T) 1 / 3 Q. Note that g u (0) is accurately 
known from the work of Gori-Giorgi and Perdew p2[ . 
The high and low frequency limits of 3fJe^4(w) are both 
obtained from the third moment sum rule. In particu- 
lar, the vanishing of ?R.eA(0) follows from the fact that 
2 joo ) . g e q ua j ( mmus ) the first moment of 

the current-current response function, which, by gauge 
invariance and the continuity equation, coincides with 
the third moment of the density-density response func- 
tion, i.e., — A(po). 

The u> 3 behavior of 5mi(w) at low frequency is easily 
obtained from the approximate zero-temperature formula 
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which is exact in the limits of high density and high fre- 
quency. Here — ^p-, and Xo-o-'(9j w ) are ^ ne s P m 
density response functions of the homogeneous liquid. 

We have evaluated x<y<y' m the generalized random 
phase approximation 

Xaa'il^) = X0a(<l,u)5<r<r> - Vf[l - G a<jl (q)] , (15) 

where xo<r(<7, w) is the Lindhard function and G,t,t' (q) are 
local field corrections p3fl . At typical metallic densities 
we multiply 3toj4(oj) by an empirical factor 



l + 



o(r a ,C) + 



(16) 



designed to satisfy the sum rule dieA(0) = with- 
out altering the high-frequency behavior. Notice that 
o(r a ,C) — > 1 f° r r s — * 0. The results evaluated with this 
procedure are shown in Figs. 1 and 2. 

Finally, we briefly remark on the calculation of the 
regular part B^ a ,(uj) of /" ccrcr ,. A simple approxima- 
tion strategy is as follows. We rewrite the 2x2 ma- 
trix fxca-a- 1 m ^ ne basis of the vectors \r) oc (nf,n±) and 
\s) oc (—71 j., rii). It is immediately seen that only the ma- 
trix element f" c rr is finite in the limit k — > 0, while all 



the others are singular. This suggests that we approx- 
imate f" c rr ~ B" r and completely ignore the contribu- 
tion from B" a ,(uj) in all other matrix elements, which are 
dominated by the singular contribution of A(u>). A little 
thought shows that this approximation is equivalent to 
setting 



(17) 



where the scalar function B a {uf) can be extracted from 
a calculation of the xc kernels in the density channel. 
Such a calculation has been carried out in p4] (for the 
paramagnetic state) and the connection between B a (uj) 
and the f£ e (w) of that paper is B a {uj) = ^ f ^ n \ y . 

This completes the construction of the input necessary 
to the evaluation of Eq. ([l0|). It is hoped that our ex- 
pression for the spin-current dependent xc field ( |Io|) will 
open the way to novel applications of CDFT to the cal- 
culation of spin excitations in spin-polarized systems. 
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